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Recent experiments at Fermilab and CERN have observed a strong asymmetry be- 
^ . tween the hadroproduction cross sections of leading D mesons, containing projectile 

valence quarks, and nonleading charmed mesons, without projectile valence quarks. 
The observed correlations of the vr^A^ D^X cross section with the projectile 
charge violates the usual assumption that heavy quark jet fragmentation factorizes. 
We examine the asymmetry between leading and nonleading charm production as 
a function of x j and pj^ assuming a two-component model combining leading- twist 
fusion subprocesses and charm production from intrinsic heavy quark Fock states. 
We predict a sizable asymmetry at low pj, and high Xf from coalescence of the 
charm quarks with the comoving spectator quarks of the projectile. An intrinsic cc 
production cross section of 0.5 /xb is sufficient to explain both the magnitude and 
kinematic dependence of the asymmetry. In contrast, the charm jet hadronization 
mechanisms contained in PYTHIA predict a sizeable leading charm asymmetry even 
at low xp. The two-component model is extended to predict the asymmetry in B 
meson production in proton-proton and pion-proton interactions. 
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I. Introduction 



In leading-twist QCD, the factorization theorem |^ predicts that the fragmentation 
functions Dh/c{z, Q) are independent of the quantum numbers of both the projectile and 
target. However, strong flavor correlations between the produced particle and the projec- 
tile have been reported in charm production p, ^. For example, in ir^iud) interactions 
with hadrons or nuclei, the D~{cd) Xf distribution is consistently harder than the D^{cd) 
distribution. The D~ and D^{cu) are referred to as "leading" charmed mesons while the 
D+ and D (cm) are "nonleading" . This leading behavior thus suggests that hadronization 
at large xp involves the coalescence of the produced charm or anticharm quarks with 
the spectator quarks of the projectile, just as in exclusive reactions. The study of this 
phenomena thus can provide new insights into the coherent mechanisms controlling the 
formation of hadrons in QCD. 

Charmed hadron distributions are often parameterized as oc (l — Xf)^, where nnonieading 
is observed to be larger than nieading [0? H]- However, this kind of parameterization is rather 
insensitive to the details of the mechanisms responsible for leading charm effects. A more 
sensitive observable, the asymmetry between leading and nonleading charm, has been 
used in the recent analyses of the WA82 and E769 collaborations. The asymmetry, 
defined as 

^ cr(leading) — a (nonleading) 
cr (leading) + a (nonleading) 

does not require the assumption of the functional form of the cross section. Both ex- 
periments find that the measured asymmetry A{xf), integrated over pt, increases from 
~ for X/ near zero to ~ 0.5 around Xf = 0.65 [||, ||. However, the asymmetry Ai^p"^), 
integrated over all Xf, is found to be consistent with zero in the range < < 10 GeV^ 
1^. These facts are consistent if the leading charm asymmetry is localized at large xp, 
involving only a small fraction of the total cross section. 

The experiment WA82 measures charmed hadron production by 340 GeV tt~ beams 
on W/Si and W/Cu targets P]. Experiment E769 has measured production in Be, Al, 
Cu, and W targets [0 with a mixed beam (75% vr" and 25% 7r+) at 250 GeV. The D+ is 
leading for a 7r+ projectile. In each experiment only mesons are used in the asymmetry 
analysis to avoid ambiguities in the assignment of leading and nonleading charm to 
mesons. The D^, a leading charm state when directly produced by vr" beams, can also 
be produced by D~^* and D^* decays; e.g. 50% of nonleading D+*'s decay to 

Perturbative QCD at leading order predicts that c and c quarks are produced with 
identical distributions. Next-to-leading order calculations do give rise to a small charge 
asymmetry (~ 10% for X/ ~ 0.8) between c and c production due to qg and qq interference 

^. However, this charge asymmetry should result in an increase of D^, D production 
over and at high Xf, not a separation between D~ , and D~^, D . 

How can one explain the origin of leading charm asymmetry within the context of 
QCD? It is clear that the produced charm (or anticharm) quark must combine with a 
projectile valence quark. Ordinary jet fragmentation {e.g. Peterson fragmentation pO| ) 
cannot produce a leading particle asymmetry since it is independent of the initial state 
and thus the projectile valence quarks. This is an essential property of leading-twist 
factorization. However, one expects on physical grounds that a charm quark produced by 
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fusion may coalesce with a comoving spectator valence quark [11, 12, 13, T^\. For example, 
in QED, leptons of opposite charge moving with similar velocities can be captured into 
neutral atoms |T5|. Since the capture is significant only at small relative rapidity, Ay, the 
effect on the total rate is higher twist. 

In leading-twist QCD heavy quarks are produced by the fusion subprocesses gg QQ 
and qq —>■ QQ. The heavy Q ot Q normally fragments independently; however, there is a 
finite probability that it will combine with a spectator valence quark in the final state to 
produce a leading hadron. Coalescence is expected to dominate when the valence quark 
and the produced heavy quark have the same velocity. The coalescence amplitude should 
be largest at small relative rapidity since the invariant mass of the Qq system is minimal 
and the binding amplitude of the heavy meson wavefunction is maximal. This picture 
of coalescence is also consistent with "heavy quark symmetry" [^, |17|. A similar final- 



state coalescence mechanism is contained in PYTHIA, a Monte Carlo program based on 
the Lund string fragmentation model [|1^]. Its string mechanism produces some charmed 
hadrons with a substantially larger longitudinal momentum than the charmed quarks 
originally produced by the fusion processes. At large Xf and low invariant string masses, 
the produced D~ or inherits all the remaining projectile momentum while D^, D 
production is forbidden. However, PYTHIA substantially overestimates the observed 
asymmetry A{xf), particularly at low xj. It also results in A{p'^) ~ 0.3 for < pf, < 10 
GeV^, overestimating the effect seen in the ^/-integrated data. 

The difficulty that PYTHIA has in reproducing the shape of the leading charm asym- 
metry seems to be a characteristic feature of models based on final-state coalescence. If 
we parameterize the fusion-produced charm distribution as dNc/dx oc (1 — x)", then the 
leading charmed hadron spectrum falls more slowly due to the extra momentum supplied 
by the projectile quark dNo/dxf oc [1 — Xf{l — 5)]" where 6 ~ mT-frnT^j and mx is 
the transverse mass. However, the asymmetries generated by this parameterization rise 
faster than the observation at low Xf that A{xf) is negligible when Xf < 0.25 |^. 
Thus a coalescence model strong enough to reproduce the hard spectrum of the leading 
charmed hadrons starting from the relatively soft fusion-produced charmed quarks tends 
to produce an excessivly large leading charm asymmetry over all phase space. It would 
be illuminating if the (true) rapidity distributions of the leading and nonleading charmed 
hadrons were measured to see whether or not their rapidity distributions track those of 
the underlying fusion-produced charmed quarks. 

In the above picture of leading charm hadropro duct ion, it is implicitly assumed that 
coalescence is strictly a final-state phenomenon. In fact, the coalescence of the charm 
quark and a projectile valence quark may also occur in the initial state. For example, 
the 7r~ can fiuctuate into a \udcc) Fock state. The most important fiuctuations occur at 
minimum invariant mass Ai where all the partons have approximately the same velocity. 
Characteristically, most of the momentum is carried by the heavy quark constituents of 
these Fock states. As viewed from the target rest frame, the intrinsic charm configurations 
can have very long lifetimes, of order r = 2Piab/ M.'^ where Piah is the projectile momentum. 
Intrinsic charm hadroproduction occurs dominantly when the spectator quarks interact 
strongly in the target explaining why large Xj charm production on nuclear targets 
is observed to have a strong nuclear dependence, similar to that of the inelastic hadron- 
nucleus cross section. Since the charm and valence quarks have the same rapidity in 
an intrinsic charm Fock state, it is easy for them to coalesce into charmed hadrons and 
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produce leading particle correlations at large xj where this mechanism can dominate the 
production rate. This is the basic underlying assumption of the intrinsic charm model 

The leading charm asymmetry must be a higher-twist effect or it would violate pQCD 
factorization. Final-state coalescence is higher twist since only a small fraction of the 
fusion-produced heavy quarks will combine with the valence quarks. Intrinsic heavy 
quark production is also higher twist-because the virtual configurations in the projectile 
wavefunction must be resolved during their limited lifetime. The cross section decreases 
with extra powers of l/mq relative to leading-twist fusion. From a general quantum- 
mechanical standpoint, both types of higher-twist mechanisms, coalescence of fusion- 
produced charm in the final state and coalesence of the intrinsic charm configurations in 
the initial state, must occur in QCD at some level. 

In this paper, we shall calculate the asymmetry within a two-component model: par- 
ton fusion with coalescence, and intrinsic charm with valence-quark recombination pT|. 



II. Leading-Twist Production 

The inclusive cross section for a single charmed hadron as a function of Xf, Xf 
[2771x1 \/s) sinhy, and in leading twist QCD has the factorized form 



^ / Hab{Xa,Xb)^ ^^/''^^^^ dz3dy2dp^dxf , (2) 
Z J H/i Z3 



where Ei is the energy of the charmed quark, y2 is the rapidity of the charmed antiquark, 
Hab{xa,Xh) is the convolution of the differential cross section with the parton distribution 
functions, and = 1.5 GeV. We shall use the lowest-order parton fusion calculation. 
A K factor of ~ 2 — 3 is included in the normalization of the fusion cross section. For 
consistency with the LO approximation, we use two current leading order sets of parton 



distribution functions, GRV LO and Duke-Owens 1.1 |2^ for the proton, and their 
pion counterparts, GRV LO ||2^ and Owens set 1 [p5[| . 

The fragmentation function, Dn/dzs), describes the hadronization of the charmed 
quark where = x^/xc is the fraction of the charm momentum carried by the charmed 
hadron, assuming it is coUinear with the charmed quark. We have studied two different 
fragmentation functions, a delta function, ^(2:3 — 1), and the Peterson function extracted 
from e~^e~ data [|I^]. We have shown previously that the Peterson function predicts a 
softer Xf distribution than observed in hadropro duct ion, even at moderate Xf. The delta 
function model assumes that the charmed quark coalesces with a low-x quark spectator 
from the sea (or a low momentum secondary quark produced in the collision) so that 
the charmed quark retains its momentum and velocity Either choice of fragmenta- 
tion function is independent of the initial state and does not produce flavor correlations 
between the projectile valence quarks and the final-state hadrons. 

In Fig. 1 we show the xj distributions calculated for (a) tt~p and (b) vr+p interactions 
at 250 GeV and (c) 7i~p interactions at 340 GeV using each set of structure functions with 
both choices of charm quark fragmentation function. The solid and dot-dashed curves 
give the calculated distributions using delta function and Peterson function fragmentation 
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with GRV LO. The dashed and dotted curves illustrate the same calculations with DO 1.1. 
For charm production at these energies, DO 1.1 gives a somewhat larger cross section. As 
expected, delta function fragmentation results in harder distributions than those predicted 
by Peterson fragmentation for Xf > 0.2. Note however that the conventional fusion model, 
even with delta function fragmentation cannot account for the shape of the leading D 
production cross section from WA82 [|]. (Since the normalization of the data has not 
been fixed, we normalize it to our calculated cross section.) The differences in the fusion 
cross sections from the GRV and DO distributions are only apparent at low Xf where 
gg cc dominates since the gluon distributions are uncertain. At large Xf, where qq — >• cc 
is more important, little difference can be observed in the xj distributions because the 
valence distributions are relatively well measured. In the following we will use the GRV 
LO distributions only. 

For valence quark coalescence to be effective as a leading charm production mecha- 
nism, either the c or c must be comoving with the projectile valence quarks since capture 
into a bound state wavefunction favors constituents with similar rapidities. Is this pos- 
sible in the fusion model? Most of the produced charmed quarks produced in the fusion 
reaction presumbably hadronize into charmed mesons or baryons independent of the pro- 
jectile identity. To first approximation, the spectator quarks have the same rapidity as the 
projectile itself since they are bound state components of the projectile wavefunction. The 
produced heavy quarks tend to have low rapidity compared to the pion valence quarks. 
In the Lund string fragmentation model, as contained in PYTHIA [|18|, a charmed quark 
is always found at the endpoint of a string. This string pulls the charmed quark in the 
direction of the other string endpoint which is usually a beam remnant. When the two 
string endpoints are moving in the same general direction, the charmed hadron can then 
be produced at larger longitudinal momentum that then original charmed quark. In the 
extreme case where the string invariant mass is too small to allow the production of several 
particles, the string scenario reduces to a coalescence one, with the two string endpoints 
determining the flavor content of the produced hadron [26]. In Fig. 2(a) we compare the 
pion valence and sea quark rapidity distributions with the fusion-produced charmed quark 
rapidity distribution at 340 GeV. The solid and dotted curves are the pion valence and sea 
quark distributions while the dashed curve is the calculated charmed quark rapidity dis- 
tribution. At this energy, the kinematic limit for charm production is reached at y = 2.8. 
The spectator valence quarks tend to have larger rapidity than the charmed quarks. In the 
region where the distributions overlap, the sea quark rapidity density is also important. 
Thus we shall assume in this paper that charmed hadrons created from the leading twist 
fusion subprocesses arise dominantly either from independent fragmentation, coalescence 
with the projectile sea components, or coalescence with comoving secondary partons pro- 
duced in the collision. Therefore we will also assume that the fusion mechanism produces 
a negligible asymmetry between and D~. We model the coalescence process with 
delta function fragmentation. 



III. Intrinsic Heavy Quark Production 

The fluctuation of a vr" into a \udcc) Fock state produces a leading particle asymmetry 
through recombination of the intrinsic cc pair with the comoving valence quarks. The 
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charmed quarks in the Fock state may be freed through soft interactions of the hght 



valence quarks with the target [|19[. The probabihty distribution corresponding to an 



n-particle Fock state (integrated over k±i) is assumed to have the form 

- iV„a,(m,s)-^ , (3) 



dxi... dxn " (ml - YJl:=i(jni / Xi)y 

where = \J {k\^ + mf is the average transverse mass, {k^) is proportional to the square 
of the quark mass, and is the normalization. Here we only consider the four-particle 
Fock stateQ. We have assumed the effective values = 0.45 GeV for the valence quarks 
and mc = 1.8 GeV for the charm quarks. 

The X distribution of intrinsic c quarks in a pion is 

^ ^ dxidx2dxc- — , (4) 



dxc J dxi . . .dx 



c 



where Xi and X2 represent the valence quark momentum fractions. The intrinsic charm 
quarks hadronize both through fragmentation as in the parton fusion model and through 
coalesence with valence quarks. Fragmentation leads to distributions of the form 



dP f dPic 

~ j dXcdZ3DD/c{z3)-J^S{xD - ZsXc) , (5) 



dx 



D 



where Xc is the fraction of the projectile momentum carried by the c quark and the 
definition of Z3 is identical to that in Eq. (2). The fragmentation of an intrinsic charm 
state itself does not produce a leading particle effect, since and production are 
equally likely. 

The coalescence of one or both of the intrinsic charm quarks with comoving spectator 
valence quarks from an intrinsic charm state naturally produces leading charmed hadrons 
at large xj. (By duality, one can think of the leading charmed mesons as preexisting in 
the intrinsic heavy quark Fock state.) The |7r^) = \udcc) state may coalesce into \uc) 
and \dc), and D~ , automatically producing leading charm since the charge conjugates 
contain no valence quarks. The leading D~ distribution is calculated from 

^"^'■^ dxc^^6{xD- - Xd- Xc) , (6) 



dxD- J dx 

where the delta function defines the sum of the c and valence d momentum 

fractions. In Fig. 2(b) we show the rapidity distributions for a light quark in a four particle 
intrinsic charm state (solid curve) and a charm quark (dashed curve). The D~ distribution 
from Eq. (6) is also shown (dot-dashed curve). The velocity (rapidity) of a heavy quark 
should remain unchanged by hadronization, up to order Aqcd/'^q- Thus there should be 
no acceleration of the heavy quark. However, since momentum is conserved, the hadron 
produced by the coalescence of equal rapidity partons will have the combined momentum 
of the heavy quark and valence spectator quark: xd- = Xg + Xc- If coalescence dominantly 



"'^Introducing additional light quarks or gluons reduces both the probability that the pion will fluctuate 
into this configuration and the probability that the cc will recombine with a valence quark to produce 
leading charm due to the presence of sea quarks. 



6 



occurs when the combining partons have the lowest invariant mass, as in this model, then 
Xq/Xc ~ mTq/fTtT-- contrast, charmed hadrons produced by jet fragmentation have 
less momentum than the parent charmed quark. Thus it is natural that leading charmed 
hadrons created by the coalescence of a produced quark with a valence quark will have 
a harder longitudinal momentum distribution than nonleading charmed hadrons. This 
effect can perhaps also explain the pattern of leading strange-charmed baryons observed 
at large x/ in the CERN WA62 hyperon beam experiment 



As in Ref. [^, we assume that the intrinsic charm model produces nonleading charm 
states by fragmentation of the intrinsic charm states and leading charm by both frag- 
mentation and valence quark coalescence. Valence quark coalescence gives rise to the 
difference between the leading and nonleading distributions. The total charm distribu- 
tion is then the sum of the parton fusion and intrinsic charm components, i.e. da/dxf = 
dapi/dxf + daic/dxf. 

Previously we obtained the normalization of the intrinsic-charm component by assum- 
ing that the ratio of the intrinsic charm cross section to the total charm cross section is 
identical to the ratio of the "diffractive" to the total J/ip production cross section mea- 
sured by NA3, giving (Tic/o"*§*^' = 0.18 for pion-induced production ||28[. The constant 
A^4 was then fixed from this ratio under the assumption that the probability distribution 



and the differential cross section are identical ||2T[. However, because of the uncertainties 
in relating cc distributions to quarkonium production, it is difficult to use the J/ip cross 
section to obtain the absolute normalization of the intrinsic charm probability. 

We set A^4 from the probability Pic, with Pic = 0.31% based upon calculations by Hoff- 



mann and Moore compared to EMC muoproduction data |29| and assume the same 
Pic for pion and proton projectiles. The intrinsic charm cross section will be proportional 
to the total inelastic cross section evaluated at s = (1 — Xf)s since a soft interaction in 



the target breaks the coherence of the Fock state and brings the particles on shell |T9|, ^ 
The cross section thus includes a resolving factor /i^/4m^, so that 



,2 



where /i^ is a hadron scale parameter. Fixing /i^ ~ 0.2 GeV^ from the NA3 ratio of 
0.18 at 200 GeV [28|, we obtain o"ic ~ 1 yub for a single charmed hadron with x/ > |. 



(The intrinsic cc pair production cross section is 0.5 /ib.) The nonleading and leading 
distributions are then 



dai 



NL 



dxf 
dah 



dx 



f 




(8) 
(9) 



where ^ is a parameter determining the relative importance of fragmentation and valence 
quark coalescence^ Following Ref. ||2l]], we have chosen ^ = 1/2 in this analysis, except 



where otherwise noted. 



^Our value of o-jc does not contradict the results of the E653 collaboration from 800 GeV pSi interac- 
tions H, adis{D+)/a{D+) < 1.8%. At 800 GeV, using cr^'^ in Eq. (7) and including an A° ''^ dependence 
for protons |||, we find cJi^{D) / (t{D) ~ 1.1%. 

•^We have not made any distinction between the relative rates of charged and neutral D production 
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The nuclear dependence, A, of the fusion model is A°' where a ~ 1 [Q. Since intrin- 
sic charm production tends to occur on the nuclear surface [|l9l, its target dependence is 
assumed to be where (3 = 0.77 for pion-induced reactions p8|. The model thus pre- 



dicts that leading charm effects are diminished in nuclear targets because of the relative 
suppression of the intrinsic charm component by A^~^. This prediction is consistent with 
the A dependence seen by the NA3 |28| and E789 |34| collaborations. We average our 
calculated asymmetry over the nuclear targets to compare with the data. 



IV. Predictions of the Two-Component Model 

The nonleading and leading Xf distributions from WA82 (circles) and the combined 
data from E769 (stars) are shown in Fig. 3(a) and (b) respectively. Because the 
data is given in arbitrary units, we give a common normalization to the calculations and 
data. Figure 3(a) shows the nonleading charm distributions for 7r~p interactions at 340 
GeV. The solid and dashed curves correspond to fusion calculations using delta function 
and Peterson function fragmentation respectively with the intrinsic charm contribution 
of Eq. (8). The Peterson function result lies well below the data, even at moderate x/. 
Figure 3(b) gives our results for leading charm from 7r~p interactions at the same energy, 
including valence quark coalescence as well as fragmentation of the intrinsic charm quarks 
as in Eq. (9). 

Figure 3(c) compares our results with the measured asymmetry at 250 (circles) and 340 
GeV (stars) as well as the combined asymmetry from both experiments (squares) |Q, |^. 
The steepest increase of A{xf) arises from Peterson function fragmentation (dashed curve) 
due to its much softer nonleading distribution at large Xf. The solid curve shows the 
delta function result. Because of the different nuclear target dependence for the intrinsic 
charm contribution, the leading charm effect diminishes in large nuclear targets. The 
dot-dashed curve shows delta function fragmentation without the average over nuclear 
target to show that the relative target dependencies of fusion and intrinsic charm do not 
strongly affect A{xf). All the calculations based on the intrinsic charm model reproduce 
the general trend of the data. Increasing the coalescence contribution in Eq. (9) by 
increasing ^ would increase the calculated asymmetry. The dotted curves in Fig. 3(b) and 
(c) show the effect of changing ^ from 1/2 to 9/10 for a proton target. The leading D xj 
distribution is marginally changed, but the asymmetry is noticeably increased, bringing 
it into better agreement with the high Xf data. This suggests that high statistics single 
D distributions combined with the asymmetry data can fix relative normalization of the 
coalescence contribution. The leading D~ rapidity distribution from the two-component 
model is shown in Fig. 2(c) for both ^ = 1/2 and 9/10. Note that the total charm rapidity 
distribution is broader than the fusion result in Fig. 2(a), especially for y > 1. 

It should be emphasized that leading-twist QCD, including leading order corrections, 
does not produce an asymmetry between leading and nonleading charm p, Both 

but have rather assumed the same proportions for intrinsic charm and fusion. One could introduce 
an additional parameter in an attempt to model production differences i.e. between leading D~ and 
nonleading D , both containing c quarks. Some differences in charged and neutral D production naturally 
occur since D*^s are commonly produced first but do not decay uniformly since D*+ and both decay 
to -D°. See Ref. ||l^ for some discussion on relative D~ , D° production by fusion due to target effects. 
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the intrinsic charm model and PYTHIA produce an asymmetry with Xf. In contrast to 
the PYTHIA predictions given in 0, ^, the intrinsic charm picture produces a slightly 
negative asymmetry for Xf < 0.4. This is due to the difference between coalescence and 
fragmentation at low Xf. The D meson Xf distribution resulting from fragmentation, 
calculated using Eq. (5), peaks at {xf) ~ 0.25, while the coalescence of the valence quark 



with intrinsic charm quarks, Eq. (6), produces leading charm with {xf) ~ 0.5 |21|. When 



Xf > 0.4, valence quark coalescence dominates over all other contributions, accounting for 
the rise in the calculated asymmetry. Without this intrinsic charm coalescence mechanism, 
i.e. ^ = 0, there would be no asymmetry in this model. The asymmetry produced by 
intrinsic charm, shown in Fig. 3(c), corresponds to the single charm production cross 
section cTjc ~ 1 yub. At higher energies, as the fusion cross section increases, the relative 
intrinsic charm contribution to the charm production cross section will diminish since cxJJJ, 
is nearly constant with energy. However, this will not affect the asymmetry at large xj 
until the fusion cross section dominates the intrinsic charm contribution. 

The effects of intrinsic charm and coalescence occur dominantly at low pt where the 
valence and charm quarks are aligned. In order to illustrate this effect, we note that the 
Pt dependence of these models is approximately 

dNpi_ ^ ag(p^ + mg) ^^^^ 

dpr [Pt + ^c]^ ' 

dNic ai{pl + ml) 

dpi + • ^ > 

The higher-power falloff in the pt dependence of intrinsic charm reflects its higher-twist 
nature. The higher power of appears in Eq. (11) because the gluon vertices with both 
the pion valence quarks and the intrinsic cc pair must be included in the calculation of 
the intrinsic charm amplitude, leading to in the cross section. In the parton fusion 
model, the structure functions include the coupling of gluons and quarks to the incident 
hadrons, thus the fusion rate is proportional to a^. 

The E769 collaboration has also investigated the p^ dependence of the asymmetry 
p. Figure 4 shows our calculated p\ dependence for several regions oi xj. Figure 4(c) 
compares the combined p^ distribution from E769 in the range 0.1 < Xf < 0.7 
1^ with the model. The sohd and dashed curves indicate delta function and Peterson 
function fragmentation respectively. We have assumed that fragmentation only affects 
the longitudinal momentum. The calculations are weighted by the percentage of the 
probability of intrinsic charm production from fragmentation and coalesence in each Xf 
region. Nonleading and leading charm calculations are shown in Figs. 4(a) and (b). The 
solid and dashed curves show our delta and Peterson function results for < xj < 0.4. 
Peterson function fragmentation, with its softer xj distribution, gives a somewhat larger 
^/-integrated cross section when < Xf < 0.4 than the delta function. Conversely, when 
0.4 < Xf < 1, the harder distribution of delta function fragmentation (dot-dashed curves) 
gives it a larger ^/-integrated cross section than the Peterson function (dotted curves), 
both for nonleading and leading charm. The difference between the delta and Peterson 
function predictions at low p^ in Fig. 4(b) is due to the relative strength of valence quark 
coalescence in the forward region. As p'^ increases, the difference caused by the higher 
power falloff in p^ is reduced primarily because of the higher power of as in Eq. (11). 
This difference will manifest itself most apparently in the asymmetry. 
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In Fig. 5(a), we compare the ^/-integrated (0.1 < Xf < 0.7) asymmetry of E769 [^j 
with our results. There is no asymmetry for > because the integrated probabihties 
for leading and nonleading charm are equal. In the region covered by the data, there 
is a slight asymmetry, shown by the solid line for delta function fragmentation and the 
dashed line for the Peterson function. Both curves are multiplied by a factor of 50 to be 
visible. The choice of the delta function leads to a negative asymmetry due to its smaller 
integrated fragmentation component from Eq. (9). The Peterson function tends to make 
the intrinsic charm fragmentation component narrower, so that most of its contribution is 
contained in the region 0.1 < < 0.7, resulting in a positive Aip^). For < x/ < 0.4, the 
calculated Aij)^) is negligibly small and negative and is not shown. When 0.4 < Xf < 1, 
the resulting A{p'^) is larger and positive, as shown in Fig. 5(b). The solid curve shows 
delta function fragmentation, the dashed, Peterson function. Both decrease strongly with 
p^, due in part to the different powers of the coupling constant in the distributions. 
Although the delta function is more physical in low p\ hadroproduction, the Peterson 
function should take over at sufficiently large p\ since it describes jet fragmentation in 
e^e~ annihilation A nonzero A{p%) is more likely to be found in the forward Xf data 
at low py. 

Recently, the Fermilab E791 collaboration has taken high precision charm produc- 
tion data with a vr" beam at 500 GeV on carbon and platinum targets in the range 
—0.1 < X/ < 0.8 ||3^. We have calculated the Xf and p^ distributions and the corre- 
sponding asymmetries for a proton target for this experiment using GRV LO structure 
functions and delta function fragmentation. The fusion cross section increases by a factor 
of two between E769 and E791, decreasing the relative strength of the intrinsic charm 
contribution by ~ 50 % in the higher energy experiment. This reduces the asymmetry 
by~ 10%atx/ ~ 0.8 for E791. The pi dependence of the asymmetry is again the 
strongest for the forward xj region, 0.4 < xj < 0.8. (Below xj ~ 0.4 A{xf) is again 
slightly negative; the predicted asymmetry becomes positive for x/ > 0.4.) We predict 
little change in A{p'^) between E769 and E791. The fusion cross section will continue 
to increase with incident energy while A{xf) will correspondingly decrease, albeit more 
slowly than the intrinsic charm contribution to the total cross section. 



V. Predictions for Leading Beauty Hadroproduction 

We also use the two-component model to predict B meson distributions. The parton 
distribution functions are more stable for bb production, thus our calculations with the 
GRV LO and DO 1.1 sets are very similar. We only show delta function fragmentation 
since there is less distinction between the Peterson and delta function predictions for b 
quarks. One surprising prediction of our model is that a larger fraction of the inclusive 
6-quark cross section is produced from intrinsic beauty compared to the fraction of open 
charm produced by intrinsic charm at fixed energy. Of course, at fixed values of r = 
rngg/v^, the fusion cross section decreases by ~ {mc/mb)'^{as{mg^) / as{mcc)Y between 
charm and beauty production. The normalization of the intrinsic heavy quark cross 
section contains factors of mQC^si'm'QQy so that the probability for producing intrinsic 
beauty relative to charm only decreases as {mc/frih)'^ rather than {mc/rfiijY as may be 
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expected from Eq. (3). An additional decrease of {fnc/rfib)'^ is associated with the resolving 
factor, yU^/4m^. We assume = 4.6 GeV and take fi"^ to be independent of both the 
projectile and the final state. Thus, the intrinsic beauty cross section is predicted to be 



on the order of 2-3 nb, similar to the size of the fusion cross section at the energies we 
investigate. The nonleading and leading B distributions are similar to those given in Eqs. 
(8) and (9). The fusion cross section is expected to increase rapidly with energy but tTib is 
proportional to the inelastic cross section. Therefore although o"ib /ct*^*^-! is relatively large 
for our calculated distributions, it will decrease at higher energies. 

No data on B distributions or their associated asymmetries have yet been published. 
We have calculated B production in 7r~p interactions at 250 and 340 GeV. Our predic- 
tions of the leading and nonleading distributions and the asymmetry between them is 
shown in Fig. 6. The solid lines in Figs. 6(a) and (c) illustrate the fusion mechanism with 
delta function fragmentation. The dashed curves show nonleading B^{hu) production 
with intrinsic h quark fragmentation as in Eq. (8) and the dot-dashed curves show leading 
S~(6u)'s as in Eq. (9). The difference between the fusion and the nonleading B distribu- 
tion is larger than for charm since Uic and Uib are coupled to the inelastic cross section. 
As the fusion cross section increases, cTib/cr*^*^^ is reduced, seen by comparing Figs. 6(a) 
and (c). In Figs. 6(b) and (d) we show the corresponding asymmetry calculations. The 
asymmetry is similar for charm and beauty because the shape of the intrinsic distribution 
is only weakly dependent on the heavy-quark mass. 

Since fixed-target experiments can be done at Fermilab with an 800 GeV proton 
beam, we include pp predictions for B production in Figs. 6(e) and (f) using the same 
notation as above. An intrinsic heavy quark state in the proton has at least a five quark 
configuration, thus a leading 5+ or a Z^^ produced by valence quark coalescence has 
{xf) ~ 0.4 while intrinsic heavy quark fragmentation gives {xf) ~ 0.2, both smaller than 
in a pion projectile. The intrinsic heavy quark cross section with proton projectiles is 
somewhat larger since cr™ is 40% larger than aj^p. The Xf distribution from fusion is 
narrower in proton production, as shown in Fig. 6(e). A faster increase in the asymmetry 
is expected (see Fig. 6(f)) since (x/) is reduced for the five-quark Fock state. 

One interesting test of the extension of this model to B production is the shape of 
the leading and nonleading distributions. Since (Jib is comparable to the bb production 
cross section by fusion, a parameterization of the xj distribution as (1 — xf)'^ should give 
a harder distribution and smaller n than expected from fusion production, especially for 
pp production, measurable at current energies. As the incident energy increases, n should 
increase as the fusion contribution becomes dominant. The difference between leading 
and nonleading B distributions may become important for GP- violation studies where an 
understanding of the symmetries of the B hadroproduction cross section is crucial. This 
is particularly true for kaon beams in order to produce BgS, with a large Lorentz gamma 
factor. 




(12) 
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VI. Conclusions 



We have shown that the intrinsic charm model produces an asymmetry between non- 
leading and leading charmed hadrons as a function of xj and p'^. The asymmetry, A, is 
predicted to be largest at low p"^ and to be an increasing function of Xf. The model also 
accounts for the shapes of the charmed hadron distributions in Xf. It would be interesting 
to look for an asymmetry between Ac and Ac production both in tt~p and pp interactions. 
Since both Ac and Ac contain a pion valence quark, the asymmetry should be small. How- 
ever, the asymmetry could be large in pp interactions where two of the proton valence 
quarks can coalesce with a charm quark to produce leading Ac while the Ac should be 
centrally produced. We also note that charmed-strange mesons are produced at large xj 



by hyperon beams a leading particle effect of the type studied here. 

The asymmetry between D~ and production in pp collisions should also be checked. 
If the two-component model is correct, the asymmetry should begin to increase at lower 
Xf than the calculated asymmetry with a pion projectile but should have approximately 
the same shape. We also predict that the B and D asymmetries should be similar. An 
additional check on our model comes from the shape of the B distributions at large Xf, 
especially with a proton beam where the intrinsic beauty contributions should produce a 
broader distribution than expected from leading twist fusion subprocesses. 
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Figure Captions 



Figure 1. The Xf distributions for (a) ir^p and (b) n^p interactions at 250 GeV and (c) 
7r~p interactions at 340 GeV calculated from parton fusion. The solid and dot-dashed 
curves show calculations using GRV LO structure functions with delta function and Pe- 
terson function fragmentation. The dashed and dotted curves illustrate the same results 
using DO 1.1 structure functions. The leading D~ data from the WA82 collaboration [^] 
are shown with the fusion calculation in (c). 



Figure 2. The rapidity density of valence (solid curve) and sea quarks (dotted) in the 
GRV LO pion distributions |Q are compared with the rapidity density of the produced 
charm quarks (dashed) in (a). The valence (solid) and charm (dashed) rapidity distribu- 
tions in an pion fluctuation into an intrinsic charm state are shown in (b), along with the 
resulting leading D distribution from Eq. (6). In (c) we show our two-component model 
calculation including fusion with the leading charm distribution from Eq. (9). The model 
calculations with intrinsic charm have been converted to rapidity. All the calculations are 
for a 340 GeV vr" beam. 



Figure 3. Our results for (a) nonleading charm and (b) leading charm distributions in 
7r~p interactions at 340 GeV and (c) the asymmetry compared with the WA82 (cir- 
cles) and E769 |^ (stars) data. The combined asymmetry from both experiments is 
also shown (squares) The calculations are with GRV LO distributions using delta 
function (solid) and Peterson function (dashed) fragmentation with the intrinsic charm 
contributions to nonleading, Eq. (8), and leading, Eq. (9), charm production. The dotted 
curve in (b) shows the leading D distribution with ^ = 9/10 in Eq. (9). In (c), the dashed 
curve is calculated with the Peterson function and the solid curve with delta function 
fragmentation. Both are averaged over nuclear target. The dot dashed curve uses delta 
function fragmentation and a proton target. The dotted curve shows the leading contri- 
bution calculated with ^ = 9/10 in Eq. (9) for a proton target. 



Figure 4. Nonleading and leading charm distributions are shown in (a) and (b) for 250 
GeV TT~p interactions. The solid and dashed curves are calculated using delta function 
and Peterson function fragmentation in the < Xf < 0.4 interval while the dot-dashed 
and dotted curves are calculated for 0.4 < Xf < 1. The combined data from E769 

is compared to our calculation integrated over 0.1 < Xf < 0.7 in (c). The solid and 
dashed curves use delta function and Peterson function fragmentation respectively. 



Figure 5. The x/-integrated (0.1 < Xf < 0.7) asymmetry of E769 is shown in (a). Our 
results, using delta (solid) and Peterson (dashed) fragmentation functions, are multiplied 
by a factor of 50. The asymmetry in the forward region, 0.4 < Xf < 1, is shown in (b). 
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The solid and dashed curves are again our results for the delta and Peterson fragmenta- 
tion functions. 



Figure 6. Our predictions of leading and nonleading B distributions are shown in (a), 
(c), and (e). The solid hues illustrate the fusion mechanism. The dashed curves show 
nonleading B production with intrinsic b quark fragmentation and the dot-dashed curves 
are leading B's {B^{bu) for Ti~{ud) and B^{ub) ior p{uud)). In (b), (d), and (f), we show 
the corresponding asymmetry calculations. 
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